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HIGHLIGHTS 


•  A  multi-scale  thermo  electrochemical  model  for  lithium  ion  cells  is  developed. 

•  The  test  results  of  lithium-ion  cells  are  presented. 

•  The  operating  condition  for  the  best  performance  of  the  examined  battery  cell  is  identified. 

•  The  results  of  applying  battery  cell  in  a  realistic  driving  cycle  is  presented. 

•  A  crucial  recommendations  are  suggested  for  the  battery  management  systems  of  BEV. 
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A  substantial  part  of  electric  vehicles  (EVs)  powertrain  is  the  battery  cell.  The  cells  are  usually  connected 
in  series,  and  failure  of  a  single  cell  can  deactivate  an  entire  module  in  the  battery  pack.  Hence,  un¬ 
derstanding  the  cell  behaviour  helps  to  predict  and  improve  the  battery  performance  and  leads  to  design 
a  cost  effective  thermal  management  system  for  the  battery  pack.  A  first  principle  thermo  electro¬ 
chemical  model  is  applied  to  study  the  cell  behaviour.  The  model  is  in  good  agreement  with  the 
experimental  results  and  can  predict  the  heat  generation  and  the  temperature  distribution  across  the  cell 
for  different  operating  conditions.  The  operating  temperature  effect  on  the  cell  performance  is  studied 
and  the  operating  temperature  for  the  best  performance  is  verified.  In  addition,  EV  cells  are  examined  in 
a  realistic  driving  cycle  from  the  Artemis  class.  The  study  findings  lead  to  the  proposal  of  some  crucial 
recommendation  to  design  cost  effective  thermal  management  systems  for  the  battery  pack. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Driving  cycles  are  a  standard  way  to  compare  the  performance 
of  different  vehicles.  There  are  two  classes  of  driving  cycle,  modal 
cycle  such  as  the  new  European  driving  cycle  (NEDC)  and  transient 
cycle  such  as  the  Assessment  and  Reliability  of  Transport  Emission 
Models  and  Inventory  systems  (Artemis).  The  NEDC  cycle  is  a 
compilation  of  straight  acceleration  and  constant  speed  periods 
representing  typical  driving  patterns,  and  uses  to  assess  the  CO2 
and  fuel  efficiency  performance  of  all  passenger  vehicles,  either 
internal  combustion  engines  or  electric  vehicles,  in  Europe. 
Whereas  transient  cycles  such  as  Artemis  involve  many  speed 
variations,  typical  of  on-road  driving  conditions  and  is  a  statistical 
representative  of  real  driving  cycles. 
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The  study  of  the  driving  cycle  has  applied  predominantly  at  the 
battery  pack  level  in  electric  vehicles.  Very  limited  research  has 
looked  at  the  effect  of  a  driving  cycle  at  the  cell  level.  Dubarry  et  al. 
built  a  real  life  drive  cycle  and  applied  a  simple  empirical  battery 
model  1  .  They  attempted  to  provide  a  coherent  approach  for  a 
realistic  empirical  model  to  simulate  battery  cell  performance,  and 
establish  a  scheme  to  better  understand  and  predict  battery  cell 
performance.  More  recently,  Klein  [2]  simplified  the  first  principle 
models  for  lithium-ion  cells  and  applied  driving  cycles  to  examine 
the  fidelity  of  the  simplified  model.  However,  the  thermal  issue  was 
not  subject  of  their  works.  Gerver  [3]  and  co-authors  studied  large 
scale  batteries  (20  x  30  cm  size)  used  under  high  power  conditions, 
5  C  rate(l  C  is  a  full  discharge  in  1  h,  based  upon  theoretical  ca¬ 
pacity,  so  5  C  is  a  full  discharge  in  12  min).  They  developed  a 
coupled  electrochemical  battery  model  with  a  thermal  model  to 
determine  temperature  distributions  throughout  the  battery.  They 
applied  their  model  to  optimise  the  terminal  tab  position  in  the  cell 
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Fig.  1.  The  model  sequence. 


and  some  similar  design  applications  of  the  cell.  However,  the 
model  results  were  compared  with  the  experimental  measure¬ 
ments  of  other  researchers  [4]  and  the  effect  of  driving  cycles  and 
heat  dissipation  parameters  were  not  studied.  More  recently,  Xiao 
[5]  and  Guo  [6]  developed  dynamic  models  for  a  pouch  type 
lithium-ion  battery  based  on  electrochemical  and  thermal  princi¬ 
ples  to  analyse  performances  of  a  single  cell.  However,  the  effect  of 
ambient  temperature  and  driving  cycle  were  not  studied  in  their 
work.  It  is  concluded  that  the  heat  transfer  issues  of  EV  battery  cells 
have  not  been  fully  incorporated  into  models  so  far  developed.  The 
heat  transfer  studies  are  focused  on  the  battery  pack  and  the  issue 
has  not  been  examined  at  the  fundamental  level  of  EV  batteries,  the 
cell. 

This  paper  is  a  development  of  a  previous  paper  published  by 
the  authors  [7]  in  which  the  results  of  study  the  battery  cell 
behaviour  in  a  typical  driving  cycle,  NEDC  were  presented.  The 
current  paper,  by  concentrating  on  the  heat  transfer  study  at  the 
cell  level,  contributes  to  understanding  the  cell  behaviour  under 


different  operating  conditions.  The  authors  developed  physic  based 
model  [8]  is  applied  to  verify  the  operating  temperature  for  the  best 
performance.  In  addition,  the  operating  temperature  effect  on  the 
cell  performance  is  studied  and  the  BEV  cells  are  examined  in  a 
realistic  driving  cycle  from  the  Artemis  class.  The  study  findings 
lead  to  the  proposal  of  some  crucial  recommendation  to  design  cost 
effective  thermal  management  systems  for  the  battery  pack. 


2.  The  model  definition  and  validation 

The  model  8],  which  is  applied  in  this  paper,  is  a  physics-based 
model  that  consists  of  two  sub-models,  the  ID  electrochemical  sub 
model,  and  2D  thermo-electro  sub  model,  which  are  coupled  and 
solved  concurrently.  The  equations  of  the  model  are  given  in  ap¬ 
pendix  A.  The  ID  model  can  predict  the  heat  generation  rate  (Qh) 
and  voltage  (V)  of  the  battery  cell  through  different  load  cycles.  The 
2D  model  of  the  battery  cell  accounts  for  temperature  distribution 
and  current  distribution  across  the  surface  of  the  battery  cell.  The 
battery  cell  considered  is  a  Lithium  Manganese  Oxide  (LMO)  cell  of 
the  pouch  structure,  and  the  cell  is  a  high  capacity  cell  (20  Ah) 
dedicated  for  electric  vehicle  application. 

In  terms  of  the  solution  sequence,  the  model  input  is  a  load  cycle 
that  is  representative  of  the  expected  drawing  down  of  power/ 
current  from  the  battery  cell  with  time.  The  ID  electrochemical 
model  for  the  load  cycle,  predicts  the  voltage  and  state  of  charge  for 
each  components  of  the  cell.  The  voltage  at  a  point  in  the  1 D  model, 
which  corresponds  to  the  tab  in  the  battery  cell,  and  is  mapped  to 
the  tab  position  in  the  2D  model.  In  addition,  the  associated  total 
heat  generation  rate  which  includes  reversible  (reaction)  and 
irreversible  (Ohmic  and  overpotential),  maps  to  each  computa¬ 
tional  cells  in  the  2D  model.  The  2D  model  is  a  thermo-electric 
model  that  solves  the  charge  balance  (Joule  equation)  and  energy 
equation  concurrently.  See  Fig.  1.  The  model  results  are  compared 
with  the  experimental  results.  Fig.  2  shows  the  model  prediction  for 
average  temperature  of  the  cell  surface  and  the  measurement  of 
average  temperature.  The  differential  equation  that  describes 


Fig.  2.  Comparison  between  the  predicted  voltage,  temperature  of  the  cell  and  the  test  results  throughout  the  load  cycle,  the  nominal  capacity  of  the  cell  20  Ah. 


A.  Tourani  et  al.  /  Journal  of  Power  Sources  268  (2014)  301—314 


303 


Fig.  3.  The  total  heat  sources  characteristic  obtained  by  the  model  at  various  discharge  rates  (1C,  20  Ah,  the  battery  cell  nominal  capacity). 


temperature  distribution  in  the  battery  takes  following  conserva¬ 
tion  form 


pCp— —  AV2T  +  qg  —  Qd  (1) 

qd  is  the  rate  of  heat  exchanges  between  the  battery  cell  and  sur¬ 
roundings  and  qg  is  the  heat  source  term  in  the  battery  cell. 


Qd 


^cell 

Eell 


(T-Ts) 


X  = 


EM 

E'i 


Thermal  conductivity  in  z  direction 


The  heat  capacity  for  the  equivalent  cell  is  volume  averaged  as 
follows 


pcp  -  y  (pcM/yji 


h  is  the  heat  transfer  coefficient  (W  m-2  K_1),T  is  the  cell  temper¬ 
ature  and  Ts  is  the  surrounding  temperature.  Aceii  and  Vceii  are  the 
surface  area  and  volume  of  the  battery  cell. 

Bernardi  [9]  governed  the  complete  heat  generation  term,  gg, 
and  Rao  and  Newman  [10]  showed  that  the  heat  of  mixing  in  the 
Bernardi  equation  is  negligible  when  the  open  circuit  voltage  (OCV) 
of  battery  does  not  change  rapidly.  Therefore,  on  provision  that  the 
reaction  distribution  is  uniform  across  the  porous  electrode,  and 
side  reactions  are  negligible,  they  simplified  Bernardi  equation  11  . 
Rao  equation  has  cited  frequently  in  the  literature  and  is  applied  in 
the  model  of  this  paper.  The  heat  source  term  (qg)  in  the  battery  cell 
is  defined  by  Ref.  [10] 


%  —  asj*nj  (0s  —  0e  —  Uj)  +  ^2  asjh M  +  o'eff^0sV0s 

j  j 

ke^V0eV0e  +  /<p)ffVln  qV0e 

(3) 

As  stated  in  Ref.  [8  ,  the  first  term  on  the  right  hand  side  is 
potential  deviation  in  the  control  volume  from  the  equilibrium 
potential  (irreversible  heat).  The  second  term  represents  the 
entropic  effects  (reversible  heat).  The  third  term  shows  the  Ohmic 
heat  in  the  solid  phase  and  the  last  two  terms  represents  the  Ohmic 
heats  in  the  solution  phase,  the  thermal  conductivity  of  the  battery 
is  anisotropic. 

Thermal  conductivity  in  X  and  Y  direction  (/;  denotes  the 
thickness  of  layer  i): 


3.  Heat  generation  in  the  cell 

The  model,  [8],  is  applied  to  verify  the  heat  generation  rate  at 
different  discharge  rates  for  0.1  C,  0.3  C,  1  C,  and  3  C  where  C  is  the 
nominal  capacity  of  the  cell  (20  Ah),  see  Fig.  3.  The  figure  shows 
the  cell  operating  at  high  current  rate  (3  C)  results  in  average  heat 


Discharge  curves 


Fig.  4.  Voltage  curves  for  the  battery  cell  at  0.3  C  discharge  rate  (6  A)  and  for  different 
ambient  temperatures. 
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Fig.  5.  The  Artemis  cycles. 


Table  1 

A  typical  electric  vehicle  parameters. 
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Frontal  /  area  of  car 

2.39 

m2 

Percent  grade 

1 

Vehicle  mass 

1720 

kg 

Drivetrain  efficiency  (inverter  to  motor) 

0.95 

Tyre  rolling  resistance  coefficient 

0.015 

Brake  and  steering  resistance 

0.003 

Drag  coefficient 

0.38 

Air  density 

1.2 

kg  m 

Wheel  radius 

0.3 

m 

Drag  coefficient 

0.38 

Air  density 

1.2 

kg  m 

Wheel  radius 

0.3 

m 

Battery  efficiency  (battery  to  inverter) 

0.95 

Gearbox  efficiency  (motor  to  wheels) 

0.95 

Number  of  cells 

300 

The  regenerating  coefficient 

0.25 

Battery  cell  capacity 

20 

Ah 

Battery  pack  energy 

25 

kWh 

generation  of  73.2  kW  nrT3  whilst  at  low  current  rate  (0.2  C)  the 
average  heat  generation  during  discharge  process  is  854  W  m-3. 
One  of  the  key  points  arising  from  the  results  is  that  the  heat 
source  in  an  EV  battery  cell  is  not  large  enough  to  satisfy  other 
heating  requirement  such  as  cabin  heating.  However,  the  heat 
generation  in  an  EV  battery  cell  application  can  be  managed  to 
maintain  the  cell  at  the  desired  operating  temperature  in  order  to 
improve  the  performance  of  the  cell.  Other  source  of  heating  in 
EVs  particularly  electrical  devices  such  as  an  inverter  and  tran¬ 
sistors  should  also  be  considered  in  designing  the  thermal  man¬ 
agement  system. 


4.  The  operating  temperature  effect 

The  lithium-ion  cell  (LMO)  is  examined  at  different  operating 
temperatures.  The  cell  average  temperature  is  considered  as  a 
design  parameter,  which  varies  from  -33  °C  to  60  °C.  The  battery 
cell  performance  is  expressed  by  voltage  discharge  curves  when  the 
nominal  current  (20  A)  is  drawn  from  the  battery.  This  cycle  is 
selected  as  representative  of  cell  performance  for  different  oper¬ 
ating  temperatures. 

As  it  can  be  seen  in  Fig.  4,  the  operating  temperature  has  sig¬ 
nificant  effect  on  the  voltage-capacity  characteristic  of  the  cell.  The 
temperature  dependent  parameters  are  valid  between  -7  °C  and 
60  °C;  however,  the  model  predicts  the  cell  voltage-capacity 
characteristic  at  different  operating  temperature  from  -33  °C  to 
60  °C.  It  is  observed  that  the  cell  cannot  be  operated  at  low  tem¬ 
perature  as  the  voltage  drops  dramatically.  For  instance,  at  the 
operating  temperature  of  -33  °C,  the  voltage  of  the  battery  de¬ 
creases  quickly  to  the  cut  off  voltage  (2.9  V)  and  only  5%  of  the 
capacity  of  the  battery  cell  can  be  drawn  down.  The  main  reason  for 
performance  dropping  is  the  reduced  rate  of  the  chemical  reaction 
associated  with  the  mass  transfers. 

At  high  temperature,  faster  chemical  reactions  occur  and  the  cell 
delivers  constant  voltage  throughout  the  discharge  cycle.  As  Fig.  4 
shows,  increasing  the  operating  temperature  from  20  °C  to  60  °C 
did  not  have  a  significant  effect  on  the  voltage-capacity  charac¬ 
teristic.  In  addition,  the  high  operating  temperatures  cause  thermal 
stress  within  the  cell  [12  .  Therefore,  an  upper  operating  temper¬ 
ature  limit  for  the  cell  should  be  considered.  In  order  to  optimise 
the  performance,  and  according  to  the  results  of  the  operating 
temperature  effect  study,  the  cell  requires  a  working  temperature 
range  of  20  ±  5  °C  and  the  cell  should  not  encounter  high  tem¬ 
perature  differences. 


5.  Drive  cycle  study 

The  fuel  consumption  and  pollutants  emissions  of  conventional 
internal  combustion  cars  are  assessed  using  standard  modal 
driving  cycles  such  as  the  new  European  driving  cycle  (NEDC).  The 
results  of  studying  the  cell  in  NEDC  cycle  was  presented  in  the 
previous  paper  published  by  the  authors  [7  .  Car  manufacturers 
use  realistic  driving  cycle  such  as  Artemis  cycle  to  better  under¬ 
standing  real  driving  conditions  and  evaluate  real  performances  of 
their  vehicles.  The  Artemis  cycle  consists  of  3  different  configu¬ 
rations,  plus  an  additional  variant.  These  are  the  urban  cycle,  the 
road  cycle,  the  motorway  130  km  hr1  and  the  motorway 
150  km  h-1,  see  Fig.  5. 

The  model  is  applied  to  investigate  the  battery  cell  behaviour 
associates  with  realistic  driving  cycles.  Typical  electric  vehicle  pa¬ 
rameters  are  considered  to  determine  the  power  required  from 
battery  cell.  The  effects  of  three  configurations  of  Artemis  cycles  are 
studied  separately  and  a  combined  motorway  and  urban  cycle  is 
also  studied,  and  further  the  effect  of  ambient  temperature  and 
heat  transfer  coefficient  are  discussed. 

Based  on  the  basic  formulae  for  the  physics  of  motion  and  for 
the  typical  battery  electric  vehicle  parameters,  Table  1,  the  power 
required  from  the  battery  cell  associated  with  the  Artemis  cycle 
load  is  calculated.  The  detail  of  the  calculation  is  described  in 
Appendix  B. 

The  power  required  from  the  battery  cell,  voltage  and  mean 
temperature  for  an  hour  journey  on  the  Artemis  urban  cycle  and 
half  an  hour  journey  on  the  Artemis  road  cycle  are  shown  in  Figs.  6 
and  7.  The  average  temperature  and  voltage  of  the  battery  cell  are 
obtained  at  different  ambient  temperatures  with  the  range 
from  -10  °C  to  35  °C. 

A  typical  electric  vehicle  is  considered  to  apply  25%  of  deceler¬ 
ation  power  to  charge  the  battery,  13]  and  the  results  show  that 
the  battery  cell  is  not  fully  discharged  either  at  the  end  of  an  hour 
urban  Artemis  or  half  an  hour  Artemis  road  journey.  This  finding 
can  contribute  in  estimating  BEV  mileage  and  demonstrates  the 
effectiveness  of  BEV  application  for  commuting  trips.  Since  a  report 
from  “Element  Energy”  [14]  identifies  commuting  as  accounting  for 
25%  of  total  annual  car  miles,  and  around  two  thirds  of  commuting 
trips  are  less  than  ten  miles  long.  This  means  that  the  battery  cell 
and  consequently  the  battery  electric  vehicle  can  achieve  the  ma¬ 
jority  of  commuting  trips  without  being  fully  discharged.  In  addi¬ 
tion,  electric  vehicles  can  apply  ultra-capacitors  [15]  or  range 
extender  technologies  to  increase  the  power  capability  or  extend 
the  range. 

The  results  of  Figs.  6  and  7  show  that  the  cell  temperature  is 
dominated  by  the  ambient  temperature  in  both  urban  and  road 
cycles.  The  cell  temperature  at  the  end  of  the  journey  reaches  the 
ambient  temperature  whether  the  ambient  temperature  is  less  or 
more  than  the  initial  temperature  of  the  cell,  20  °C.  This  is  assuming 
a  heat  transfer  coefficient  of  h  =  6  W  m-2  K_1  for  edges  and  large 
surfaces  of  the  cell.  This  heat  transfer  coefficient  can  be  obtained 
from  the  different  methods  of  heat  transfer,  convection,  conduction 
and  radiation.  This  finding  can  lead  to  verify  the  objective  of  the 
thermal  management;  as  it  would  be  concentrating  on  insulating 
the  cell  from  adverse  ambient  temperature. 

The  cell  behaviour  is  also  studied  for  a  half  an  hour  journey  on 
the  motorway  Artemis  cycle;  Fig.  8.  The  results  are  in  agreement 
with  the  results  of  using  the  battery  cell  in  urban  and  road  Artemis 
cycles  as  the  cell  is  not  fully  discharged  at  the  end  of  journey  and 
the  cell  temperature  is  again  dominated  by  ambient  temperatures. 
However  since  the  voltage  decreases  to  the  cut  off  voltage,  the  cell 
could  not  manage  to  last  longer  than  half  an  hour  journey  with 
Artemis  motorway  cycle.  In  addition,  the  cell  temperature  rises 
above  the  ambient  temperature  and  the  configured  thermal  mass 
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Fig.  6.  The  battery  cell  required  power,  voltage  and  mean  temperature  in  an  hour  journey  with  Artemis  urban  cycle  for  different  ambient  temperatures. 
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and  heat  dissipation  rates  could  not  balance  with  the  heat  gener¬ 
ation  rates. 

In  order  to  investigate  the  battery  cell  behaviour  in  a  realistic 
journey  between  two  cities,  a  typical  combined  journey  within 
Artemis  motorway  and  urban  cycles  is  considered.  It  is  assumed 


that  the  journey  is  combination  of  6  min  urban  driving,  maximum 
speed  50  km  hr1,  followed  by  18  min  motorway  driving,  maximum 
speed  120  km  h  \  and  6  min  urban  driving.  The  power  required, 
voltage  and  average  temperature  of  the  battery  cell  are  shown  in 

Fig.  9. 
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Fig.  7.  The  battery  cell  required  power,  voltage  and  mean  temperature  for  a  half  an  hour  journey  with  Artemis  road  cycle  for  different  ambient  temperatures. 
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Fig.  8.  The  battery  cell  required  power,  voltage  and  mean  temperature  in  half  an  hour  journey  with  Artemis  Motorway  cycle  for  different  ambient  temperatures. 
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Fig.  9.  The  battery  cell  required  power,  voltage  and  mean  temperature  in  half  an  hour  combined  cycle  (Artemis  Motorway  and  Urban  cycles)  for  different  ambient  temperatures. 
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The  results  show  the  voltage  and  average  temperature  of  the  cell 
change  in  a  way  to  reach  the  ambient  temperatures  for  the  first 
6  min  of  the  urban  cycle,  then  the  voltage  drops  and  temperature 
rises  up  gradually  through  the  motorway  cycle.  The  journey  ends 


up  with  6  min  urban  cycle  during  which  the  cell  temperature  falls 
down  to  the  ambient  temperature. 

The  study  again  shows  the  effect  of  the  ambient  temperature 
and  highlights  that  the  heat  dissipation  rate  should  be  increased  in 


cr> 

CL> 


Fig.  10.  The  battery  cell  mean  temperature  in  half  an  hour  Artemis  Road  cycle  with  constant  heat  transfer  coefficients  and  different  ambient  temperature. 
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motorway  segment  of  the  cycle.  This  finding  leads  to  designing  the 
thermal  management  of  BEV  cells  in  a  way  to  provide  variable  heat 
dissipation  rates  depending  on  the  heat  generating  rates  and  the 
ambient  temperatures. 

6.  The  heat  transfer  study 

Battery  cells  depending  on  their  position  in  the  battery  pack 
might  have  different  heat  transfer  coefficients,  and  surrounding 
conditions.  In  this  study,  different  possible  heat  dissipation  condi¬ 
tions  are  considered  for  a  particular  battery  cell,  which  in  com¬ 
parison  with  thermal  studies  of  a  battery  pack  are  less  complex,  and 
lead  to  lower  solution  times  and  cost.  Instead,  the  multi  physics 
based  model  for  different  load  conditions  predicts  the  heat  gener¬ 
ation  term.  The  heat  transfer  coefficient  is  considered  as  a  design 
parameter  with  the  range  0,  6, 16  W  m”2  K-1,  and  the  cell  behaviour 
at  different  ambient  temperature  and  heat  transfer  coefficients  is 
studied.  The  heat  transfer  coefficient  of  zero  corresponds  to  an 
adiabatic  condition  in  which  the  battery  cell  is  assumed  to  be 


Ambient  Temperature  5degC 


insulated  from  the  surrounding.  The  heat  transfer  coefficient  of 
6  W  m~2  K_1  corresponds  for  a  passive  cooling  with  a  buoyancy 
flow  around  the  battery  cell,  and  the  highest  heat  transfer  coeffi¬ 
cient  represents  active  cooling  with  a  forced  flow  around  the  bat¬ 
tery  cell.  The  Artemis  road  cycle  with  the  average  speed  of 
57  km  h-1  is  considered  as  the  load  cycle  and  the  results  are  shown 
in  Figs.  10  and  11.  The  temperature  of  surroundings  or  equally  a 
coolant  temperature  can  represent  the  cell  ambient  temperature,  if 
active  cooling  is  applied. 

Fig.  10  shows  the  results  of  using  different  heat  transfer 
coefficients  within  the  constant  ambient  temperatures.  It 
is  observed  when  the  cell  is  in  the  adiabatic  condition, 
h  =  0  W  nrr2  K~\  the  cell  temperature  eventually  increases  by 
3  °C  in  30  min  regardless  the  variation  of  the  ambient  tempera¬ 
ture.  The  cell  in  the  adiabatic  condition  with  the  considered 
thermal  inertia  would  cope  with  the  amount  of  heat  generated 
throughout  the  cell. 

As  the  cell  begins  to  transfer  heat  with  the  surrounding,  h  =  6, 
16  W  m~2  I<-1  in  Fig.  10,  the  cell  experiences  higher  temperature 


Ambient  Temperature  -lOdegC 


Fig.  11.  The  battery  cell  mean  temperature  in  30  min  Artemis  Road  cycle  with  different  heat  transfer  coefficients  and  constant  ambient  temperature. 
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gradient  than  the  adiabatic  conditions  and  this  gradient  depends 
on  the  ambient  temperature.  For  example,  the  cell  with  the  highest 
heat  transfer  coefficient,  16  W  m-2  K-1,  and  the  lowest  ambient 
temperature, -10  °C,  has  the  temperature  difference  of  30  °C 
throughout  the  half  an  hour  journey  with  Artemis  road  cycle. 

The  results  for  the  upper  limit  of  the  heat  transfer  coefficient, 
16  W  m-2  I<-1  in  Fig.  10,  shows  that  with  this  heat  transfer  co¬ 
efficient,  and  the  assumed  thermal  mass  of  the  cell,  the  cell 
temperature  is  dominated  by  the  ambient  temperatures.  On  the 
other  hand,  the  thermal  inertia  of  the  cell  is  considered  as  a 
constant  parameter  that  can  be  a  function  of  temperature  or  state 
of  charge.  Although  this  assumption  is  a  well-accepted  for  re¬ 
searches  on  the  cell  level  [16],  this  is  should  be  verified 
experimentally. 

The  studies  of  using  different  ambient  or  coolant  temperatures 
are  shown  in  Fig.  11.  The  ambient  temperature  varies  from  -10,  5, 
20, 35  °C.  For  the  ambient  temperature  of  20  °C,  using  the  high  heat 
transfer  coefficient  helps  to  reduce  the  temperature  difference 
throughout  the  cycle  and  it  appears  the  difference  drops  from  3  °C 
in  adiabatic  condition  to  1  °C  when  the  heat  transfer  coefficient  is 
16  W  m-2  K-1.  However,  it  would  not  be  a  cost  effective  approach  to 
apply  high  heat  transfer  coefficient  at  this  ambient  temperature. 
Conversely  results  have  been  shown  when  the  ambient  tempera¬ 
tures  are  -10,  5,  35  °C  using  the  high  heat  transfer  coefficient  of 
16  W  m-2  K-1  increases  the  cell  temperature  difference  from  4  °C  in 
the  adiabatic  condition  to  15  °C  when  the  ambient  temperature  is 
35  °C.  This  finding  shows  that  the  optimum  coolant  temperature 
can  be  20  °C  when  the  ambient  temperature  is  not  equal  to  this 
temperature. 

7.  Conclusion 

A  first  principle  model  is  applied  to  study  the  lithium-ion 
(LMO)  cell  behaviour  at  different  operating  conditions.  The  re¬ 
sults  show  that  below  an  operating  temperature  of  7  °C,  a  com¬ 
mon  ambient  temperature  in  a  Northern  European  winter,  the  cell 
performance  falls  off  dramatically,  whereas  there  is  little  or  no 
gain  in  the  performance  above  27  °C.  Further,  the  lithium-ion  cell 
(LMO)  is  examined  in  Artemis  realistic  driving  cycles,  and  the 
effect  of  ambient  temperature  on  the  cell  behaviour  during  the 
driving  cycles  is  studied.  The  results  show  the  BEV  cell  temper¬ 
ature  in  the  urban  cycle  is  a  function  of  ambient  temperature,  and 
the  main  thermal  concern  in  the  urban  driving  cycle  is  to  main¬ 
tain  the  cells  in  an  optimum  working  temperature  envelope 
(20  ±  5  °C). 

As  the  BEV  cell  needs  to  draw  less  than  a  1  C  discharge  rate  in 
urban  cycles,  the  cell  can  dissipate  the  heat  through  natural  con¬ 
vection.  However,  the  cell  experiences  a  rising  temperature  when 
drawing  more  than  3  C  discharge  during  high  speed  driving  or 
acceleration.  An  active  cooling  system  will  be  needed  in  this  case  to 
reduce  the  temperature.  In  addition,  the  results  highlight  the  role  of 
a  maximum  power  limit  on  the  cell  temperature,  which  can  be 
minimised  to  reduce  the  need  for  thermal  management.  Further 
analysis  is  needed  to  determine  the  optimum  number  of  cells  to 
reduce  the  power  limit,  minimise  the  cost,  weight,  size,  and  ther¬ 
mal  management  cost  of  the  battery  pack. 

The  results  of  studying  the  effect  of  heat  transfer  coefficients 
show  that  when  the  cell  is  in  an  adiabatic  condition,  the  cell 
temperature  eventually  rises  3  °C  regardless  of  the  variation  in 
the  ambient  temperature.  The  cell  temperature  at  the  end  of  the 
journeys  approaches  the  ambient  temperature  whether  the 
ambient  temperature  is  less  or  more  than  the  initial  temperature 
of  the  cell,  20  °C.  These  findings  show  the  main  thermal  issue  of 
an  urban  BEV  cell  is  insulating  the  cell  from  adverse  ambient 
temperature.  In  addition,  high  heat  transfer  coefficients 


particularly  in  a  low  temperature  surrounding  would  result  in 
cooling  of  the  cell. 

In  summary,  this  paper  presents  thermal  behaviour  study  of  a 
battery  cell,  which  is  the  preliminary  and  essential  study  to  design  a 
thermal  management  system  for  a  battery  pack,  and  helps  to  un¬ 
derstand  the  scope  of  the  thermal  problem  within  the  battery  cell. 
The  knowledge  could  apply  to  configure  the  battery  module  and 
pack  in  a  way  to  reduce  the  cost  and  weight  of  thermal  manage¬ 
ment  systems  of  BEVs.  The  cell  model  is  planned  to  contribute  in  a 
BEV  thermal  management  systems  model  as  proposed  by  Ref.  [17] 
in  a  way  that  the  cell  level  model  will  be  developed  by  an  assembly 
of  ten  cells  with  integration  of  cooling  system. 

The  findings  of  the  project  give  a  methodology  for  determining 
the  heat  generation  and  temperature  distribution  within  lithium- 
ion  cells  in  a  BEV,  and  introduce  a  new  hypothesis  to  design  a 
thermal  management  system  for  BEV  batteries.  The  thermal  man¬ 
agement  system  should  be  designed  in  a  way  to  maintain,  for  the 
cell  under  consideration,  a  varying  temperature  of  20  ±  5  °C 
whatever  the  ambient  conditions. 

It  is  recommended  that  the  study  of  heat  generation  and 
dissipation  in  EV  cells  in  battery  packs  is  continued  and  the 
performance  of  thermal  management  techniques  to  maintain  the 
cell  in  the  required  operating  range  are  compared.  For  example, 
phase  change  materials  [18]  and  heat  pipes  [19]  have  been 
introduced  recently  to  the  thermal  management  of  BEVs,  and  the 
model  can  be  modified  to  include  these  passive  cooling  methods. 
In  addition,  the  ways  to  provide  a  coolant  in  an  active  thermal 
managements  system  with  low  temperature  variation  should  be 
studied. 


Appendix  A.  The  model  electrochemical  equations 


List  of  symbols 


Cl 

Concentration  of  the 
electrolyte  (mol  m  3) 

T 

Temperature  (K) 

Cs 

Surface  concentration  in 
spherical  particle  (mol  m'3) 

F 

Faraday's  constant  (C  mol'1) 

CP 

Specific  heat  (J  kg'1  K'1) 

f 

Ionic  activity  coefficient 

jloc 

Transfer  current  density 

(A  m  2) 

4 

Transport  number 

Q 

Volumetric  heat  generation 
rate  (J  m  3  s'1) 

e 

Porosity  (volume  fraction 
of  electrolyte) 

4 

Equilibrium  potential  (V) 

Df 

Effective  diffusivity  in  the 
electrolyte  (m2  s'1) 

U-H 

U-H 

<u 

Effective  solution  conductivity 
(Q  1  m'1) 

io 

Exchange  current  density 

(A  m'2) 

A 

Thermal  conductivity 
(W  m'1  K'1) 

4ef 

Electrode  particle's 
equilibrium  potential  (V) 

P 

Density  (kg  m'3) 

Cl,  surf 

Particle  surface  concentration 
(mol  m  3) 

&eff 

Effective  matrix  conductivity 
(Q  1  m'1) 

k  o 

reaction  rate  constant 

#1 

Potential  in  the  matrix 
phase  (V) 

Cl,  max 

Maximum  surface 
concentration  (mol  m  3) 

@2 

Potential  in  the  solution 
phase  (V) 

V 

Activation  overpotential  (V) 

k 

Electronic  conductivity 
(S  m  l) 

£>i 

Electrolyte  salt 
diffusivity  (m2  s'1) 

5a 

Specific  surface  area 

(m2  m'3) 

9ln//6ln  Ci 

Activity  dependence 

R 

Gas  constant  J/(  mol  K) 

Dsn 

Negative  electrode 
intercalation  diffusivity 

(m2  s'1) 

y 

Fraction  of  lithium  in  the 
positive  electrode 

Uref, 

Open  potential 
of  the  electrode  (V) 

T 

Bruggeman  coefficient 

X 

Fraction  of  lithium  in  the 

negative  electrode 
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k,  =  0.0001c,  (  -  10.5  +  0.668  x  10~3q  +  0.494*  10~6cf 
+0.074T  -  1.78  x  10~5q  x  T  -  8.86  x  10  10cf 
xT  -  6.96  x  10-5T2  +  2.8  x  10~8q  x  T2)  (A-14) 


D,  =  10  ^  x  10 


54 

-4  43- , 

-4  w  m  V  7  —  229  -  5  x  10_3c, 


-0.22  x  10~3q 
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Ds„  =  1.4523  x  10- 13  exp 
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Appendix  B.  Calculating  the  power  required  from  a  battery 
(A- 5)  cell 


where 

V  =  -Eref(CbSUrf) 


*0 


^Cs,max 


^s,surf  )  ^s,surf 


kf  =  k i=l,2 
Df  =  Die\  i=l,2 


(A-6) 

(A-7) 

(A-8) 

(A-9) 


Three  major  forces  which  resist  a  vehicle  from  moving  are: 
rolling  resistance,  aerodynamic  drag  and  the  force  of  gravity  as  a 
vehicle  moves  up  a  hill.  These  forces  are  calculated  and  the  battery 
cell  required  power  from  the  cell  is  estimated  for  a  typical  electric 
vehicle,  Table  1.  In  addition,  25%  of  the  total  required  power  is 
considered  to  recharge  the  battery  during  deceleration  [13] 


^roll  —  Froll^H?  (B  1) 

Froii:  the  rolling  force 
/iron:  the  rolling  resistance  coefficient 
M:  the  mass  of  the  vehicle  (kg) 
g :  the  force  of  gravity  (9.8  m  s~2) 


^ref,n[  ] 


[iref,p[21] 


-0.16  +  1.32  exp(  -  3x)  +  10  exp(  -  2000x) 

(A-10) 

-23.53y5  +  52.43y4  -  44.86y3  +  18.4y2  -  3.35y 
+  3.16 

(A-ll) 


D=\/2CdAPV2  (B-2) 

D:  the  drag  force 
C+  the  Drag  coefficient 
A:  the  frontal  area  of  the  car  (m2) 
p:  the  density  of  the  air  (kg  m-3) 

V:  the  speed  of  the  vehicle  (m  s-1) 

The  force  required  to  move  uphill  is  a  function  of  the  hill  angle 


From  Ref.  [22] 


Fhin  =  (%  grade  of  hill )Mg 


(B-3) 


344.1347148  exp(-32.9633287x  +  8.31671 1484)/(1  +  749.0756003  exp(-34.7909964x  +  8.887143624)) 
-  0.8520278805X  +  0.36229929x2  +  0.2698001697 


(A-12) 


(4.131274309  exp(0.571536523y)  +  1.281681122  sin(-4.9916739y)  -  0.090453431  sin(-20.9669665y  +  12.5788250) 

-  .0313472974  sin(31.7663338y  -  22.4295664)  -  4.14532933  +  8.1471 13434y  -  26.064581y2  +  12.7660158y3 

-  0.184274863  exp(  -  (y  -  0.5169435168)/0.04628266783)2) 


(A-l  3) 
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The  power  required  to  move  a  vehicle,  Pm,  at  a  certain  speed  is 
calculates  by 

^m  =  (Froll+D  +  F_hill)V  (B-4) 

The  power  required  to  accelerate,  Pa,  from  the  speed  of  V\  to  V2 
in  the  time  interval,  t,  is  calculated  by  the  kinetic  energy  equation. 


The  total  power  is  the  summation  of  the  acceleration  and 
movement  powers. 
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